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ABSTRACT 



New stellar models which track the production and destruction of ^He (and D) have been 
evolved for a range of stellar masses (0.65 < M/Mq < 100), mctallicitics (0.01 < Z/Zq < 1) 
and initial (main sequence) ^He mass fractions (10^^ < X-i,^ms ^ 10^'^). Armed with the ^He 
yields from these stellar models we have followed the evolution of D and "^He using a variety 
of chemical evolution models with and without infall of primordial or processed material. 
Production of new ^He by the lower mass stars overwhelms any reasonable primordial 
contributions and leads to predicted abundances in the presolar nebula and/or the present 
interstellar medium in excess of the observationally inferred values. This result, which 
obtains even for zero primordial D and ^He, and was anticipated by Rood, Steigman & 
Tinsley (1976), is insensitive to the choice of chemical evolution model; it is driven by the 
large ^He yields from low mass stars. In an attempt to ameliorate this problem we have 
considered a number of non-standard models in which the yields from low mass stars have 
been modified. Although several of these non-standard models may be consistent with the 
^He data, they may be inconsistent with observations of ^^C/^^C, ^^O and, most seriously, 
the super-^He rich planetary nebulae (Rood, Bania & Wilson 1992). Even using the most 
extreme of these non-standard models (Hogan 1995), we obtain a generous upper bound to 
pre-galactic ^He: X3P < 3.2 x 10~^ which, nonetheless, leads to a stringent lower bound to 
the universal density of nucleons. 
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1. INTRODUCTION 



To test and constrain models of primordial nucleosynthesis it is necessary to confront 
the predicted abundances with observational data. Such comparisons are complicated by 
the necessity of extrapolating the abundances inferred "here and now" (solar system, ISM, 
etc.) to their "there and then" (primordial or pregalactic) universal values. This difficulty is 
somewhat ameliorated for ^He which is observed in less evolved, low metallicity extragalactic 
HII regions (e.g., Pagel et al. 1991; Skillman et al. 1994; Olive & Steigman 1995) and for 
^Li which is probed in metal poor halo stars (Spite & Spite 1982; Thorburn 1994). In 
contrast, ^He is only observed in the solar system (Black 1971, 1972; Geiss & Reeves 1972; 
Geiss 1993) and in Galactic HII regions (Bania, Rood & Wilson 1987; Balser et al. 1994). 
Until recently, D, too, had only been observed "here and now". A new era in probing 
deuterium has dawned with the identification of possible absorption of QSO light by high 
redshift D (Carswell ct al. 1994; Songaila et al. 1994; Tytler 1995). However, "wrong 
velocity" hydrogen absorbers may masquerade as deuterium (Carswell et al. 1994; Songaila 
et al. 1994; Steigman 1994) so that much more data are required before such observations 
may be used to fix the nearly primordial abundance of D. Thus, to infer the pregalactic 
abundances of D and ^He it is currently necessary to have recourse to models of galactic 
chemical evolution. Unfortunately, the uncertainties and/or lack of uniqueness of such 
models compound the observational uncertainties. 

The evolution of D is straightforward since when incorporated in a star it is burned (to 
^He) during the pre-main sequence evolution. If the "virgin" fraction of the ISM (either 
today or at the time of formation of the solar system 4.5 Gyr ago) were known, the primor- 
dial D abundance could be inferred from ISM or solar system observations. A very large 
class of chemical evolution models (Audouze & Tinsley 1974; Tosi 1988 a,b; Vangioni-Flam 
& Audouze 1988; Mattcucci & Prangois 1989; Steigman & Tosi 1992), constrained by heavy 
element abundances, abundance gradients, abundance ratios (primary vs. secondary) and 
cosmochronomctcrs, find that ~ 1/3 — 2/3 of the ISM (now, to o^j at the time of the so- 
lar system formation, Iq) has never been through stars. However, "designer models" with 
larger D destruction (Vangioni-Flam, Ohve & Prantzos 1994; Ohve et al. 1995; Scully & 
Olive 1995) do exist (but see, e.g., Edmunds 1994 and Prantzos 1995 for their consistency 
problems) . 

The galactic cvohitiou of ^He is complex. Any prestellar ^He is enhanced by the pre- 
main sequence burning of prestellar D. Thus, when a star reaches the main sequence the 
^He mass fraction exceeds that in the prestellar nebula [X^^ms = (-^3 + 3X2/2)pre; in what 
follows we write Xj for mass fractions and yi for ratios by number to hydrogen; yij = yi+yj]. 
^He survives nuclear burning in the cooler outer layers of stars, is destroyed in the hotter 
interiors and, especially in low mass stars, is synthesized via hydrogen burning. Thus, 
depending on the mass (and, to a lesser extent the metallicity) of a star, ^He is preserved, 
produced or destroyed. To understand the galactic evolution of '^He, it is necessary to first 
understand its stellar evolution (Iben 1967; Rood 1972; Dearborn, Schramm & Steigman 
1986; Vassihadis & Wood 1993). However, since some ^He survives stellar processing and. 
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prestellar D is burned to ^He, the evolution of ^He (or, D + ^He) is less rapid than that 
of D alone (Yang et al. 1984; hereafter YTSSO). Thus, observations of D and "^He (e.g.. in 
the solar system) may be used to infer an upper bound to primordial D and/or D + "^He 
(YTSSO; Walker et al. 1991; hereafter WSSOK). As a result, models which destroy D more 
efficiently are to some extent constrained by the requirement that solar system and/or ISM 
^He not be overproduced. 

Rood, Steigman & Tinsley (1976; hereafter RST) first included stellar production of ^He 
in numerical models of galactic evolution. RST found that stellar production led to a rapid 
increase in X3 from to to which seemed in conflict with the then current observations 
(only bounds) and concluded that it seemed unlikely that '^He could be used cosmologically. 
As we shall see, this problem persists today (Ohve et al. 1995). 

In a previous study (Steigman & Tosi 1992; hereafter ST92) we used the ^He survival 
from Dearborn, Schramm & Steigman (1986) in models of galactic evolution consistent with 
a large number of observational constraints (Tosi 1988a) to track the evolution of D and ^He. 
Not including ^He production in low mass stars, we found that X3P ^ X^q ~ X^q for a wide 
range of choices of initial (primordial) abundances for D and ^He. Similar results have been 
found, e.g., by Fields (1995). At the same time, the ^He HII region data (Bania, Rood & 
Wilson 1987; Rood, Bania & Wilson 1992) have been supplemented by new and important 
observations of super-^He rich material in a few planetary nebulae (Balser et al. 1994, Rood 
et al. 1995, hereinafter RBWB95). Since the gas in PNe reflects the chemical composition 
of the ejected outer layers of their central stars, where the original ^He abundance has 
certainly been modified by stellar processing, their large ^He abundances should not be 
taken as representative of the ISM abundances at the time of the formation of their central 
stars but, of the effect of stellar nucleosynthesis in low mass stars. The observed abundances 
lend support to the estimates of newly synthesized ^He (X3* ^ 0.7— 7 x 10~^) in low mass 
stars (0.6 — 2.3Mq) in the pioneering work of Iben (1967) and Rood (1972) and makes 
timely a reinvestigation of the galactic evolution of ^He (along with D) including stellar 
production of ^He (Galli et al. 1995, Tosi, Steigman & Dearborn 1995, Olive et al. 1995). 
To this aim, we have computed an extensive grid of stellar evolution models with varying 
initial abundances of ^He and of the overall metallicity Z, and we have reexamined the 
galactic evolution of D and ^He in the framework of such stellar models. 

In §2 we describe the stellar models we've evolved to follow the evolution of ^Hc and 
present our results for the survival/production of ^He as a function of stellar mass (0.65 < 
M/Mq < 100), metaUicity {Z = 0.02, 0.002, 0.0002) and initial ("main sequence" ^D + 
^He) ^He (IO^Xs^ms = 1-0, 5.0, 10, 21, 100). In §3 and §4 we use the stellar resuhs in a 
series of chemical evolution models (Tosi 1988a,b; ST92) to follow the galactic evolution of 
D and ^He. In §5 we further follow the ^He evolution modifying the stellar yields. Sections 
6 & 7 are reserved for discussions and conclusions. 
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2. STELLAR EVOLUTION OF =^He 



The evolution of stars with M = 0.65Mq to M = IOOMq and three metahcities ( "Pop I" : 
X = 0.70, Z = 0.02, "Pop I 1/2": X = 0.76, Z = 0.002, "Pop IF: X = 0.76, Z = 0.0002) 
was followed using the code described in Dearborn, Gricst & RafFelt (1991) which is derived 
from that of Eggelton (1967, 1968). OPAL opacities (Iglesias & Rogers 1990, 1992) are used 
for temperatures above 6000K and Los Alamos opacities (Heubner et al. 1977), including 
the contribution of molecules, for lower temperatures. Our code follows ^He production 
and destruction in all hydrogen burning regions and a post-processor is used to calculate 
the nucleosynthesis of D and the isotopes of Li, Be, C, N, O and several other elements. 
Calculations with this post-processor confirmed the pre-main sequence conversion of D to 
^He. The mixing length is chosen to give the correct solar radius and this model predicts a 
neutrino flux (for the chlorine experiment) of 7.9 SNU and the p-mode spectrum matches 
that of Bahcall & Ulrich (1988). For the Pop 11/2 and Pop II models, a slightly larger 
mixing length is required to fit the color-magnitude data for M30 (Bolte 1994). However, 
^He evolution is insensitive to any reasonable choice of mixing length. 

All models began in the pre-main sequence on the Hayashi track. The final evolutionary 
state depended on the mass of the star. Stars with 0.65 < M/Mq < 2.00 were evolved up to 
the helium flash. For one model (l.SM©) the evolution was followed through the helium flash 
and up the asymptotic giant branch (AGE). This additional evolution (without thermal 
pulses) changed the ^He mass fraction in the envelope by less than 1%, demonstrating that in 
the absence of thermal pulses, the envelope abundance of ^He is insensitive to the subsequent 
(post helium flash) evolution of these stars. Intermediate mass stars (2.5 < M/Mq < 8.0) 
were all evolved to the AGB, to a point where mass loss (and possibly Planetary Nebula 
formation) dominates their remaining evolution. All massive star models (M > IOM0) were 
evolved into carbon burning and are within a few thousand years of core collapse, leaving 
insufficient time for significant changes in the composition of their hydrogen envelopes. We 
have not allowed for winds and this suggests that, for the more massive stars, we have 
underestimated the '^He returned to the ISM (see, e.g.. Dearborn, Schramm & Steigman 
1986). 

Since D is completely burned to '^He in the pre-main sequence, our "main sequence" (MS) 
^He abundance reflects the prestellar contribution from D and ^He: X^^ms — Xs + 3X2/2. 
To explore the sensitivity of our ^He yields (envelope fractions Xs^env)) to X^^ms we have 
chosen five values of X^^ms from 1 x 10^^ to 1 x 10^^. The top panel of Figure 1 shows 
our Pop I yields X^j = X3{env) as a function of the initial stellar mass for different choices 
of Xs^MS- For lower mass stars production dominates and, when convection on the giant 
branch homogenizes the envelope, the ^He available to be returned to the ISM is strongly 
enhanced. Without substantial mass loss, these lower mass stars remain near the Hayashi 
track during helium burning and beyond. In the absence of a hot bottom convection zone 
stimulated by thermal pulses, there is no destruction of this "new" ^He. For more massive 
stars the lifetimes are shorter resulting in less ^He production in the temperature regimes 
where this occurs. Note that if these stars started with no initial ^He they would, in fact. 
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be net producers of ''He; in this case = X^f / X^^ms would diverge. However, for the 
range of X^^ms of interest (10^^ — 10^"^), the small production of ^He in more massive stars 
does not compensate for the destruction of initial '^He. 

A striking feature of the 15 and 2'oMq Pop I models is the dip in X3/. The high 
envelope opacity allowed these models to homogenize the envelope on the giant branch 
prior to helium burning. The blue loop that occurred during helium burning then led to a 
second epoch of ^He destruction. In contrast, the 50 and IOOMq models completed their 
helium burning before evolving to the red thus leading to the "recovery" of X3J seen in 
the top panel of Figure 1. Similarly, the lower metallicity (Z=0.002 and Z=0.0002) 25M0 
models did not evolve all the way to the red giant branch prior to helium burning and, 
therefore, these stars only experienced ^He destruction on the main sequence (see bottom 
panel of Fig.l). For Pop 11/2 and Pop H stars the lower opacity results in a more compact, 
warmer structure (similar to that of a slightly more massive Pop I star). At a fixed mass 
this leads to a shorter lifetime resulting in less production and more destruction of ^He. 
However, the lower metallicity 15 & 25M0 models do not experience a second epoch of ^He 
destruction. As a result, they show less ^He destruction than the Pop I models (i.e., no dip 
in X^j vs. X^^Ms)- By comparing the two panels of Fig.l, it is apparent that the amount 
of ^He ejected by stars depends much more on its main sequence ^He abundance than on 
the overall stellar metallicity. 

In Figure 2 we display our Z=0.02 (Pop I) results from a different perspective by showing 
the final yield {X^f) as a function of the initial abundance {X^^ms) foi' 1 ^ M/Mq < 100. 
For M = IMq, production dominates and X^f is (nearly) independent of X^^ms (except 
at the highest X^^ms)- In contrast, for M = lOOM© production is negligible and X3/ 
varies linearly with X^^ms- The transition from low to high mass models is seen for M 
= 3 and 8Mq. For 0.65 < M/Mq < 2.5, where production dominates, we find X^f ^ 
{Mq/ MY'"^ X 10~^. For more massive stars the relation between X^f and X^^ms is more 
complex. As X^^ms increases, ^He production is relatively less important resulting in the 
curvature seen in Figure 2 for M = 3 and 8 Mq. For high masses where production is 
negligible, X3/ 0.33X3,m5- 

In the absence of '^He production it is interesting to consider the '^He "survival fraction" 
^3 = Xsf/Xs^MS (Dearborn, Schramm & Steigman 1986; YTSSO; ST92). In the presence of 
^He production is less useful and, potentially misleading since the yield is not proportional 
to X3 Afs- Nonetheless, for comparison with previously published results, g^ vs. M for our 
five choices of X3 m5 is shown in Figure 3. The divergence of the curves at low M refiects 
the increasing importance of ^He production which sets in first (at the largest M) for the 
lowest values of X^^ms- Note that in the limit that X^^ms goes to zero, gs diverges. In 
Figure 4 is shown (g^), the value of g^ averaged over Tinsley's (1980) IMF as a function of 
the lower mass limit m/ of the IMF. For gas incorporated in a generation of stars, as time 
increases material is returned from lower mass stars and so Fig. 4 provides a picture of the 
time-evolution of the ^He "survival". Note that for all rrii, (g^) ^ 0.3. Depending on the 
initial abundance, production dominates, {gs) ^ 1, for mi up to 7Mq. In terms of galactic 
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evolution this implies that the ISM starts to be enriched in ''He as soon as stars less massive 
than 3 — TMq (depending on the initial '^He abundance) start to die. We note again that 
stars above ~ 25M0 are likely to experience mass loss which will return unburnt "^He to 
the ISM. For massive stars our neglect of mass loss results in an underestimate of g^,. 

The trends in our results are easy to understand. During the hydrogen burning phase 
^He is produced in the core of stars of all masses. However, with time, the core temperature 
(which depends on stellar mass) increases to values where ^He is burned. Thus, new ^He 
survives only in a radiative shell adjacent to the convective core and, from this shell, may 
later be dredged up to the surface where it will mix with the prestellar {MS) '^He. The 
competition between net production or destruction then depends not only on X^^ms but, on 
the size of the shell and the amount of new ^He there, and on the depth of the dredge-ups 
(i.e., do they bring to the surface material which is enhanced or depleted in ^He). In low 
mass stars much new ^He is produced and most of it survives to dominate over (reasonably 
small values of) X^^ms- In contrast, for massive stars little new ^He is synthesized in the 
core, destruction is efficient and the ^He preservation shell is small so that the initial ^He 
dominates. 

We note that the Eggleton type code we use permits larger time steps on the AGB 
which tends to suppress thermal pulses. However, when short timesteps are enforced, these 
pulses occur. Such pulses are plausible sites for the lithium enhancements - produced at 
the expense of ^He - observed in some S-type stars. Although it doesn't require much ^He 
depletion to yield huge lithium enhancements, since the lithium is fragile it is possible that 
the lithium enhancement attains a steady state with ^He processed to ^He. A limit on 
such ^He processing might follow from the fact that a very fragile nucleus, ^^O, survives 
the dredge up in normal carbon stars (Dearborn 1992). Since ^He is destroyed at higher 
temperatures than those for ^^0, the lack of ^^O depletion suggests that ^He is not strongly 
depleted in normal carbon stars. 

A similar argument applies to the main sequence mixing proposed to explain the low 
12q^13q ];atios observed in red giants (Dearborn & Eggleton 1976; Dearborn 1992). In red 
giants ^^O appears to be independent of the ^"''C enhancement, perhaps due to stabilization 
by a molecular weight gradient. Since it is difficult to mix with a region of depleted ^He 
without modifying ^^O, it is suggestive that ^He is not destroyed in these stars. In contrast, 
J-type carbon stars have CNO-equilibrium values for ^^C/^^C and do show ^^O depletions. 
In such stars ^He is likely destroyed throughout the envelope. If most carbon stars pass 
through such a stage (rather than these representing a separate class of stars), little "^He may 
survive. Pinsonneault (Private Communication) has noted that for solar metallicity stars, 
the open cluster M67 seems to be the dividing line between stars which do not undergo giant 
branch mixing and those which do suggesting that mixing for stars with M ^ 1.3 — l.SM© 
is inconsistent with the ^^C/^^C data. 

It is interesting to compare our results with the recent work of Vassiliadis & Wood 
(1993). We both agree that in the low mass stars there is no significant change in ^He 
between the first and second dredge up. Further, Vassiliadis & Wood calculate thermal 
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pulses and, for stars below 5 solar masses, they find no change in the '^He abundance. Our 
mixing length approximations are nearly the same and they included a wind which stripped 
the envelope on the AGB while we did not. Quantitatively, our yields and theirs are in 
excellent agreement (provided that their 6'(^He) is a mass fraction) despite the fact that 
their opacities are different (higher) than ours and they used slightly older nuclear reaction 
rates. 



3. THE GALACTIC EVOLUTION OF D 

As in ST92 we have followed the evolution of D and ^He for the two "best" chemical 
evolution models for the galactic disk identified in Tosi (1988a) along with a third model in 
which, for comparison, infall is absent. Model 1, the "best" model (Tosi 1988a; Giovagnoli 
& Tosi 1995), consistent with the major observational constraints, has an exponentially 
decreasing SFR (with a 15 Gyr e-folding time), depending on both the gas and total mass 
density currently observed in each ring, a constant (in time) , uniform (in space) infall rate 
of 0.004 MQ/kpc^/yr and uses Tinsley's (1980) IMF. Model 25, the "second best" model 
(Tosi 1988a) also has an exponentially decreasing SFR with, however, a 5 Gyr e-folding 
time, an effectively constant (e-folding time of 100 Gyr) infall rate of uniform density 
0.002 Mq/Upc^ /yr and Tinsley's (1980) IMF. Our comparison No Infall (NI) model uses 
the same IMF and SFR e-folding time as Model 1 but is normahzed so as to reproduce 
the current SFR and gas/total mass distributions with galactocentric distance. We recall, 
however, that such a model does not reproduce the major features (distribution with time 
and galactocentric distance of the chemical abundances) observed in the galactic disk. 

Normally, our models adopt 13 Gyr for the present epoch and 8.5 Gyr for the formation 
of the solar system. To explore the sensitivity to the age of the model, we have run some 
models (indicated by a subscript 10) where the present epoch is 10 Gyr and solar system 
formation is at 5.5 Gyr. 

Since infall plays an important role in Models 1 & 25, it is necessary to specify the 
chemical composition of the infalhng gas. We have considered models with primordial 
infall {Zinf — 0) for which X2inf — and X^i^f — as well as models with partially 
processed infalling material. In the latter cases we have adopted Zinf = 0.2Zq, a value low 
enough to preserve the infall dilution efficiency and keep the model predictions in agreement 
with the observational constraints in the solar neighbourhood and in the whole disk (Tosi 
1988b, Matteucci & Frangois 1989). When Zi^f = 0.2Zq, the infall abundance of D is 
certainly lower than primordial due to stellar processing, whereas the ^He abundance in 
principle can be either lower or higher depending on the mass of the stars contributing to 
the enrichment of the accreted gas. We have thus considered models with X2inf = 0.7X2p 
or X2inf = O.8X2P and X^inf = O.8X3P or X^inf = X^p or X^inf = I.2X3P. ^From a 
series of model checks we find that larger deviations of the infall abundances of D from the 
primordial value would be rather improbable with an overall metallicity Zi^f — 0.2Zq. 
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To explore the sensitivity of our results to the primordial values we have run each model 
for a range of choices of X2P and X^p (0 < lO^Xap < 9.0; < lO^Xsp < 4.0). In Table 1 
we list the set of models explored. 

The evolution of deuterium is straightforward since any D incorporated in stars is de- 
stroyed. Thus, the D survival factor, /2 = X2/X2P, is identical to the fraction of the ISM 
that has never been through stars. /2, then, is independent of X2P and purely reflects the 
chemical evolution model. In Figure 5 we show /2 evaluated at the solar ring (R — Skpc) 
as a function of time. For each model (1 &: 25) we show the differences between primordial 
infall {Zinf = 0, X2inf = -^2p) and non-primordial infall {Zi^f = 0.2Zq, X2P = O.8X2P) as 
well as the corresponding no-infall (NI) model. We also compare in Figure 5 models whose 
present age is 13 Gyr with 10 Gyr models. The models with non-primordial infall have less 
D refuehng of the ISM and show a steeper decrease with time of /2 than models with pri- 
mordial infall. The NI model destroys D more slowly at first but, without replenishment of 
ISM D via infall, eventually has the largest D depletion. The short lifetime models (10 Gyr) 
have less time to destroy D but, to reproduce present observations have higher initial mass 
and SFR than the 13 Gyr models. The net result of this balancing act is that the 10 Gyr 
(infall) models have higher f2Q but lower /20; the NI models have larger /2 correlating with 
lower age. For Models 1 & 25 the entire ranges are: 0.49 < < 0.73; 0.43 < /20 < 0.62; 
for the NI models: 0.70 < f2Q < 0.78; 0.30 < /20 < 0.36. Thus, for our range of models 
D astration is modest, typically by a factor of 1.3 to 2.0 at ^0 and a factor of 1.6 to 3.3 
at to. Our Models 1 & 25 have very little evolution in /2 from Iq to to (0.9 < /2O//20 < 
1.2 for to — 13 Gyr; 1.3 < /2O//20 < 1-5 for to — 10 Gyr). This is entirely consistent with 
the solar system and ISM data: X20/X2Q — 1.6 ± 0.6 (e.g., Steigman & Tosi 1995). 

In Figure 6 we show the time evolution of X2 in the solar ring for Models 1, 25 & 

NI for to = 13, 10 Gyr. Also shown in Figure 6 are the 2a ranges of the solar system 
(Geiss 1993; Steigman & Tosi 1995) and ISM (Linsky et al. 1992; Steigman & Tosi 1995) 
D abundances. For Models 1 & 25 the ISM constraint dominates, limiting primordial D to 
X2P < 6 X 10~^ (for Yp < 0.25, y2p ^ 4 x 10~^). Even for the NI model we find a restrictive 
bound X2P < 9 X 10"^ {y2P < 6 x 10"^). 

In Figure 7 is shown the predicted radial distribution of D/H at present. The positive 
gradient is a natural consequence of larger D destruction in regions with larger SFR. The 
"data" are from Wannier (1980) but have been divided by a factor of 100 to facilitate 
comparison; while they appear to reflect the expected radial distribution, they are some 
two orders of magnitude too large, possibly due to the chemical fractionation of deuterated 
molecules. 

With the adoption of a chemical evolution model (e.g., 1, 25 or NI; to = 13 or 10 Gyr; 
X2inf/ X2P = 1.0 or 0.8) the evolution of deuterium (with time and location in the Galaxy) 
is completely determined. Comparison with solar system and ISM observations then fixes 
(bounds) the primordial D abundance. For our "normal" models (not the NI models) we 
find 2.7 < 10^X2P < 6.6 (for Xp > 0.75, 1.8 < 10^y2P ^ 4.4) which, for standard big bang 
nucleosynthesis (Walker et al. 1991; Thomas et al. 1995) bounds the nucleon-to-photon 
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ratio r) : 4.1 < 7710 ^ 7.1 (7710 = lO^^r] = lO^'^riN/n^). Our (artificial) NI models prefer 
somewhat higher primordial D: 4.7 < lO^Xap < 8.8 (3.2 < 10''y2P ^ 5.9; 3.5 < r]io < 5.1). 
The modest astration we find is similar to many previous results (Audouze & Tinsley 1974; 
ST92; Galli et al. 1995; Fields 1995). Nonetheless, it may be possible to construct chemical 
evolution models with more destruction of D although consistency with all the observational 
constraints is quite improbable (Vangioni-Flam & Audouze 1988; Vangioni-Flam, Olive & 
Prantzos 1994; Olive ct al. 1995; Olive & Scully 1995). An important test of all models is 
the evolution of ^He (RST). 



4. THE GALACTIC EVOLUTION OF ^He: STANDARD MODELS 

The evolution of ^Hc is much more complex than that of D since ^He may be destroyed, 
preserved and produced in differing proportions in stars of differing masses. Here, we 
have used the stellar results described in §2 in concert with the large number of models 
summarized in Table 1. The key difference with our earlier work (ST92), which utilized 
similar models, is our allowance here for the production of new ^He synthesized in lower 
mass stars. Following the results of §2, in running the models, at each time step the adopted 
stellar yields arc those corresponding to the X^f of the dying stars. As anticipated by RST, 
the effect of '^He production is large, dominating the evolution of ^He. 

In Figure 8 the evolution of '"^Hc is shown for a wide variety of choices of X2P, X^p, X2inf, 
^3inf for Model 1 and a present age of 13 Gyr. The vertical bar at 8.5 Gyr corresponds 
to the 2(7 range of abundances derived in the solar system (Geiss 1993). Given the lack 
of ^He determinations for the local ISM, the vertical bar at 13 Gyr corresponds to the 2a 
range of abundances derived by RBWB95 for HII regions between 6.4 and 10.3 kpc from 
the galactic center (but excluding W3 which definitely lies outside the average distribution 
and is the worst case for pressure broadening corrections). In the upper panel of Figure 8, 
the problem of excessive ^He production identified by RST is clear. Using the stellar yields 
from §2, it may be seen that there are no choices of X2P and X^p consistent with the range 
of ■'He abundances inferred from the solar system and/or ISM data. After the first few 
Gyr of evolution the contribution from newly synthesized ''He overwhelms any primordial 
D+'^He and, even when X2P = X^p = 0, the "'He abundances after 8.5 (13) Gyr are in 
excess of the solar system and ISM upper bounds. 

The trends displayed in the upper panel of Figure 8 are easy to understand. For fixed 
choices of X^p, X^lt) increases with increasing X2P since any prestellar D is burned to 
^He, enhancing the main sequence abundance of ^He. Models which begin with higher X^p 
always have higher ^He abundances although, with time, the differences are reduced by the 
emerging dominance of the newly synthesized ^He. 

The upper four curves in the bottom panel of Figure 8 demonstrate the effect of the 
different choices for the infall abundances. However, these differences are small compared 
to the overall enhancement by stellar produced ^He. As in the upper panel, even in the 
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absence of any primordial D or "^He, the predicted ^He abundances at the time or formation 
of the solar system and/or at present are in excess of the observational upper bounds. 

Thus, our results confirm - with a vengeance - the RST identified problem of ^He 
overproduction when the contribution of new ^He from low mass stars is included in models 
of galactic chemical evolution (see also Galli et al. 1995). This conclusion is not modified 
when we used Model 25, the NI model or, for any of these models with a 10 Gyr disk age. 
We have, therefore, considered how we might have to modify the contribution of new ^He 
from the low mass stars in order to reconcile ^He evolution with the observational data. 



5. THE GALACTIC EVOLUTION OF =^He: NON-STANDARD MODELS 

As outhned in Table 1, in addition to our standard models (with the yields from §2), we 
have considered several non-standard models by modifying the ^He yields from low mass 
stars. As may be seen from the bottom panel of Figure 8 and in both panels of Figure 9, 
there are a variety of possible solutions to the problem of "^He overproduction. The rapid 
increase in ^He is most curtailed in models VI which follow the suggestion of Hogan (1995) 
that in stars less than 2.5 Mq, ^He is destroyed on the giant branch before it can be returned 
to the ISM and the final envelopes of these stars therefore only contain an abundance 
X^f corresponding to the equilibrium value ^Hc/H=lxlO~^. However, Hogan's (1995) 
suggestion is in confiict with the normal ^"^C/^^C ratios observed in stars more massive 
than 1.3-1.5 Mq (Pinsonneault, Private Communication), as well as the observations of 
lithium (more fragile than ^He) in some of them. Even worse, the observations of excess 
=^He (Xg of order 10"=^) in three planetary nebulae (Rood, Bania & Wilson 1992; RBWB95) 
appear to confirm that stars around 1.5 M© are efficient ^He producers, in confiict with 
Hogan's proposal. Charbonncl (1995) has recently argued that the deep convective mixing 
responsible for the "^He destruction in Hogan's suggestion takes place in stars experiencing 
the helium fiash and gives an upper mass limit of 2 Mq to the phenomenon. This would 
reconcile the deep mixing ^He destruction with the high PNe abundances if the initial 
mass of the PNe progenitor was at least 2 Mq. However, the three well studied PNe 
have presumably originated from stars of initial mass 1.5 Mq (Stanghellini, 1995 private 
communication), thus lying in the range of stellar masses which should deplete and not 
enhance ^He in the deep mixing hypothesis. 

To try to retain some of the benefit manifest in models VI, we have allowed ^He to be 
reduced to its equifibrium value in lower mass stars (< 1 Mq in models V, as suggested by 
Pinsonneault, 1995, private communication; and < 1.6 Mq in models VII). However, since 
such low mass stars have long lifetimes, the reduction in ^He is only effective during recent 
epochs; solar system "^He is still overproduced in both cases V and VII (see bottom panel 
of Fig. 8 where the two curves are so close that only one is shown). Even assuming a total 
^He destruction in stars below 1 Mq (i.e. ^3^=0, case IV) does not solve this inconsistency, 
since it provides results indistinguishable from case V. 
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Any other case, with upper mass cutoff for Hogan's ^He destruction intermediate be- 
tween 1 and 2.5 would either be inconsistent with the large abundance observed in PNe (if 
the cutoff is larger than 1.6 Mq) or inconsistent with the "low" solar system abundances 
(if it is smaller than 2 Mq). We have therefore tested some models assuming that ^He is 
reduced to its equilibrium value in stars with mass between an arbitrary lower cutoff and 
2.5 Mq. If the lower mass cutoff is around 1.3 Mq (case VIII), smaller mass stars still 
produce large ^He consistent with that observed in PNe and the ISM abundances predicted 
at the various epochs are consistent with the corresponding observed values (see Figs 8 and 
9). 

In the models labelled II and III we have - arbitrarily - ignored new '^He production 
and set ga = 1 for M below 2 Mq (case II) or, between 1 and 2 Mq (case III; see, e.g., 
Wasserburg, Boothroyd & Sackmann 1995). As may be seen in Figure 8, models III are 
consistent with solar system and ISM data provided that the primordial abundances of D 
and/or ^He are not too large. 

Although some of the non-standard models II- VIII may avoid overproduction of presolar 
and/or ISM ^He the spatial distribution of "^He observed in galactic HII regions (Bania, 
Rood & Wilson 1987; Rood, Bania & Wilson 1992; Balser et al. 1994; RBWB95) provides 
an important constraint on all such models. The data are puzzling (see, e.g., Olive et al. 
1995), exhibiting no very well defined trend of ^He/H with galactocentric distance R (see 
Figure 10). Indeed, in contrast to the theoretical expectation that ^He/H should decrease 
with R (where the SFR is highest - in the inner galaxy - ^He/H should also be highest), 
the data hint at the opposite trend (see Figure 10) typical of elements like H and D which 
are destroyed and not produced by stellar nucleosythesis. The only models which avoid 
a decreasing ^He vs. R relation are those with case VI yields, but cases II and, perhaps, 
III and VIII may also be consistent with the observed radial distribution, once all the 
uncertainties are taken into account. 

We note that, despite the apparent conflicts of Hogan's (1995) suggestion with the 
12q^13q^ lithium and planetary nebulae data, comparison between case VI and the obser- 
vational data still leads to a significant upper bound to primordial '^He. For case VI we find 
that Xsq/ATsp = 1.4 so that for X3Q < 4.5 x 10"^ (Geiss 1993; ST95), X3P < 3.2 x 10"^ 
(yap < 1.4 X 10-5). 



6. DISCUSSION 

The evolution of deuterium depends solely on the chemical evolution model, and the 
D survival fraction, /2(t), is independent of the primordial abundance and of any stellar 
uncertainties since D is fully destroyed during the pre-main sequence evolution. For our 
"best" models (1 & 25; Tosi 1988a,b), D is destroyed by only a modest factor (1.4 - 2.0) 
by the time of formation of the solar system and, by a slightly larger factor (1.6 - 2.4) up 
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to the present epoch. Although the artifical NI models permit somewhat more destruction 
(by a factor of 2.8 - 3.3) by the present epoch, they actually destroy less D (1.3 - 1.4) up to 
the time of the formation of the solar system. The solar system and ISM data (Geiss 1993; 
Linsky et al. 1992) may be used along with models 1 & 25 (for a present age of either 10 
or 13 Gyr) to bound (at the 95% CL) the primordial D abundance: 2.7 < IO^Xsp < 6.6 
(1.8 < 10^y2P < 4.4). Comparing to the predictions of BBN (e.g., WSSOK) permits us to 
bound the universal ratio of nucleons to photons: 4.1 < rjio < 7.1. For our more extreme 
NI models we find: 4.7 < IO^Xsp < 8.8 (3.2 < IOV2P < 5.9) and 3.5 < r]io < 5.1. 

In contrast to D, the evolution of ^He is much more sensitive to the details of the 
chemical evolution model and, especially, to the physics of stellar structure and evolution. 
A key component in our analysis here has been the computation of a new and extensive 
grid of stellar models covering a wide range of masses (0.65 < M/M© < 100), heavy 
element abundances (0.01 < Z/Z© < 1) and initial (main sequence) ^He abundances (10~^ < 
^3,MS ^ 10^^). We have then followed the ^He evolution for models 1, 25 & NI with several 
choices for the present age of the disk (10 & 13 Gyr) and for the chemical composition of the 
infaUing gas (0.7 < Xsinz/Xsp < 1.0; 0.8 < X^inf/X^P < 1.2). As anticipated by RST, after 
the first few Gyr of evolution, stellar production of ^He dominates (see, e.g.. Figure 8). This 
contribution from newly synthesized ^He quickly overwhelms the primordial contribution 
and, even in the absence of any primordial ^He (and D), we predict solar system and ISM 
^Hc abundances in excess of those inferred from the observational data (see also Galli et 
al. 1995, Tosi ct al. 1995). Clearly, there is a problem and, until this confiict is resolved, 
^He cannot serve as a probe of BBN (RST). The problem may lie with the observational 
data and/or its interpretation. Or, it could be that our models - or some of the ingredients 
therein - are the culprits. 

The solar system ^He data (meteorites, lunar soil and rocks, solar wind) have recently 
been reanalyzed by Geiss (1993) and by Copi, Schramm & Turner (1995). Although the 
size of the error estimates has increased compared, e.g., to those used in ST92, the central 
values remain unchanged. Consistency with our model predictions would require that the 
presolar ^He abundance has been underestimated by more than a factor of 1.5-2.0. 

The ISM data are more problematic. If, indeed, stellar production of new ^He is oc- 
curring as indicated by the planetary nebulae observations (Rood, Bania & Wilson 1992; 
RBWB95), then X3 should be higher where there is more stellar processing - in the in- 
ner galaxy. As may be seen in Figure 10, there is no observational evidence for such a 
trend. Indeed, the highest ^He abundances are derived from data for the HII regions in 
the outer galaxy. The possibility that the ^He abundances inferred from HII region radio 
observations are not reliable indicators of the current ISM values has also been invoked. 
Olive ct al. (1995) have shown that unresolved structures in the nebulae may lead to an 
underestimate of their actual ^He content but, only by a few tenths. On the other hand, 
the occurrence of a strong, ^He depleted, Wolf-Rayet wind can reduce the ^He abundance 
inside the HII region, although we presume that the external layers of many HII regions 
should still be uncontaminated and that a negative internal gradient should appear in the 
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abundances derived from central to outer parts within the nebulae observed with sufficient 
spatial resolution. Besides, if an HII region exhibits depleted '^Hc because of the pollution 
from the central Wolf-Rayet star, it should also show enhanced "^He, which does not seem to 
be the case for the RBWB95 sample. Ohve et al. (1995) have tried to understand the ^He 
distribution as a reflection of the mass of the HII regions (more destruction of ^He in the 
more massive HII regions). This radical explanation makes little sense to us (and, to them 
as well) since it would require that HII regions are more efficient processors of interstellar 
material than the galaxy as a whole (since more than half the gas in these HII regions would 
have had to be cycled through the massive stars of the individual regions) . Since the HII 
region data are hard to acquire and difficult to analyze (see, e.g., RBWB95), work on both 
fronts is important and, beyond the scope of our analysis here. So, here we have adopted 
the extant data, assumed that the inferred abundances are correct and investigated the 
implications for our models. 

As may be seen in the Figures 8 & 9, the contribution from newly synthesized ^He is 
large compared to the differences among the various model assumptions regarding the age 
of the disk, the infall abundances, the primordial abundances, the specific models (1, 25, 
NI). Therefore, the prime suspect must be our estimates of the ^He yields from low and 
intermediate mass stars (Rood 1972 ; Galli et al. 1994; Hogan 1995; Wasscrburg, Boothroyd 
& Sackmann 1995). To explore this avenue we have considered a series of non-standard 
alternatives to our standard models. Indeed, several suggestions have been published for 
physical mechanisms to suppress the overproduction of stellar ^He. They are related to 
low-energy resonances or to deep convective mixing (e.g. Galli et al. 1994; Hogan 1995; 
Wasserburg, Boothroyd & Sackmann 1995), but none of them seems fully consistent with 
all the available data. In our cases II, HI, VI and VIII the suppression of stellar produced 
^He may be sufficient to flatten the X3 vs. t relation enough so that consistency with 
the local data may be found provided that the primordial abundances are small enough. 
However, it must be emphasized that if the detections of excess ^He in planetary nebulae 
reported by Rood, Bania & Wilson (1992) and RBWB95 are confirmed, some of these 
alternatives (II and VI) are excluded. By comparing the model predictions with all the 
available constraints, we find that only case HI is sufficiently compatible with all the data, 
despite its shallow negative gradient in the X3 vs. R relation. In this case we may use 
the solar system data to bound the primordial abundance of ^He from above and, the 
nucleon-to-photon ratio from below. 



7. CONCLUSIONS 

We have tracked the evolution of the abundances of D and ^He in a variety of chemical 
evolution models which incorporate the results of a newly computed grid of stellar structure 
and evolution models. We have confirmed that D is fully destroyed during the pre-main 
sequence evolution and, therefore, its galactic evolution is simple. For the range of our 
models (see §3) we find only modest destruction of D and, using the solar system and 
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ISM data in concert with our "best" models, we bound the primordial (pre-galactic disk) 
abundance, 

2.7 < 10%p < 6.6, 1.8 < 10^y2P < 4.4. (1) 

For consistency with the predictions of standard BBN (WSSOK), we require that the uni- 
versal ratio of nucleons (baryons)-to-photons lie in the range, 

4.1 < 7710 < 7.1. (2) 

In terms of the baryon density parameter, fls, {^b^^Iq = O.OlSr/io; Hq = SOhso km/s/Mpc), 

0.06 < Qfih^o < 0.10. (3) 

For our NI models the X2P estimate is somewhat higher and the bounds on the baryon 
density slightly lower. 

Our stellar models, which account for the evolution of ^He (destruction, survival, pro- 
duction) reveal that for low mass stars the production of newly synthesized ^He is very 
important for the evolution of galactic ^He. Indeed, for all of our standard models (§4) pro- 
duction of new ^He is so dominant that even in the absence of any pre-galactic D and/or 
^He, ^He is overproduced compared to the solar system and ISM data. Until this conflict is 
resolved it is difficult to see how ^He can be used as a probe of BBN (RST). Setting aside 
the possibility that the problem lies with the observational data, we have explored a series 
of non-standard models (§5) in which some or all of the newly sysnthesized ^He is assumed 
to be destroyed before being returned to the ISM. We find that if this suppression of "^He 
occurs only in stars < 1 Mq, the depletion is "too little, too late" to resolve the discrepancy 
with the solar system data. In contrast, '^He destruction in stars with masses in the range 
1-2 M0 (Wasserburg, Boothroyd & Sackmann 1995) can reconcile our models with the data 
provided that the initial ^He abundance is not too large. In this case we recover a lower 
bound to the baryon density which is consistent with that derived from the D evolution 
data (see eqs. 1-3). It remains to be seen whether this "fix" will be permitted by the PNe 
and or ^^C/^^C data. 
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FIGURE CAPTIONS 



Figure 1. Top panel: final envelope abundance of ^He as a function of the stellar 
initial mass for five different main sequence abundances of ''He and Z=0.02. Bottom panel: 
final envelope abundance of "'He as a function of the stellar initial mass for three different 
metallicities and X^^^ms = 2.1 x 10~^. 

Figure 2. Final envelope abundance of ■'He for Z=0.02 in stars of initial mass M/Mq—1, 
3, 8 and 100 as a function of the main sequence abundance. 

Figure 3. Survival fraction g-^ of ^He as a function of the stellar initial mass for Z=0.02 
and five different main sequence abundances of "'He. 

Figure 4. Average convolved with Tinsley's (1980) IMF as a function of the lower 
mass limit nii on the IMF integral. 

Figure 5. Deuterium survival factor /2 as a function of time predicted by different 
models of the solar ring (R= 8 kpc). Models shown in the left panel assume 13 Gyr for the 
current age of the galactic disk; those in the right panel assume 10 Gyr. In both panels, 
dotted curves correspond to no infall models, solid curves to models 1 (upper curve) and 25 
(lower curve) with X2inf — O.8X2P, dashed curves to models 1 (upper curve) and 25 (lower 
curve) with X2inf = ^2P- 

Figure 6. Evolution in the solar ring of the D abundance. Vertical bars give the 2a 
range for the abundances derived from solar system and local ISM observations (see text). 
Left panels correspond to a disk age of 13 Gyr, right panels to 10 Gyr. Models in the 
top panels assume the maximum value of primordial D consistent with the data; models in 
the bottom panels the corresponding minimum value. The dash-dotted lines correspond to 
"standard" models with primordial infall (1-B-Ia and lio-B-Ia), the other symbols are as 
in Fig. 5 

Figure 7. Current radial distribution of the D abundance resulting from the same models 
as Fig.5 (same symbols as in that figure). In each pair of curves of the same line-type, the 
steeper one corresponds to model 1 and the flatter one to model 25. Data points are from 
Wannier (1980) but arbitrarily divided by 100. 

Figure 8. Evolution of the ^He abundance in the solar ring as predicted by models 1. 
The vertical bars show the 2a range for the solar system and ISM abundances (see text for 
details). In the top panel all models adopt the standard stellar yields (see Sec. 2 and Table 
1) and assume primordial infall. They are in order of decreasing ^He and D abundances 
(i.e. from top to bottom): 1-R-la, 1-J-la, 1-W-Ia, 1-1-Ia, 1-X-la, 1-H-la, 1-M-la. The 
models in the lower panel all start with the same primordial abundances of D and ^He but 
have different adopted stellar yields and infall metallicity. From top to bottom they are: 
l-W-(Ic, la, lb. Id, Vb, Vlllb, Illb, VIb). 

Figure 9. Evolution of the ^He abundance in the solar ring as predicted by different 
chemical evolution models (1, 25, NI) using the same stellar yields and infall metallicity. 
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The left-hand panel corresponds to a disk age of 13 Gyr and the right-hand panel to 10 
Gyr. The lower and upper solid lines are models 1-T-IIIb and 25-T-IIIb respectively (see 
Table 1); the lower and upper dashed lines are models 1-H-VIIIb and 25-H-VIIb; the dotted 
line is for NI-Z-III. 

Figure 10. Predicted current radial distribution of the "^He/H number ratio. Data points 
and error bars are from the HII region analysis by RBWB95 and the specific models are 
labelled. 
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Table 1. Models and Initial Abundances. 



Model^ lO^Xsp lO^Xsp Model^ lO^Xsp lO^Xsp 



Models with standard ^He evolution and primordial infall 



1-M-Ia 


0.0 


0.0 


1-K-Ia 


4.0 


1.0 


1-A-Ia 


2.5 


2.0 


1-0-Ia 


4.5 


0.0 


1-H-Ia 


3.0 


0.0 


1-C-Ia 


5.0 


2.0 


1-I-Ia 


3.0 


2.0 


1-X-Ia 


6.0 


0.0 


1-B-Ia 


3.0 


2.5 


1-W-Ia 


6.0 


2.0 


1-J-Ia 


3.0 


4.0 


1-R-Ia 


6.0 


4.0 


1-L-Ia 


3.85 


0.0 


lio-B-Ia 


3.0 


2.5 


Models 


with standard 


3He 


evolution and non-primordial infall 


1-K-Ib 


4.0 


1.0 


1-W-Ib 


6.0 


2.0 


1-K-Ic 


4.0 


1.0 


1-W-Ic 


6.0 


2.0 


1-K-Id 


4.0 


1.0 


1-W-Id 


6.0 


2.0 


1-C-Ib 


5.0 


2.0 


1-W-Ie 


6.0 


2.0 


Models 


with non-standard 


^He evolution and primordial infall 


1-T-IIIa 


6.0 


1.5 


25-T-IIIa 


6.0 


1.5 


lio-T-IIIa 


6.0 


1.5 


25io-T-IIIa 


6.0 


1.5 


Models with non-standard ^He evolution and non- 


-primordial infall 


1-H-VIIIb 


3.0 


0.0 


1-Q-IIIb 


6.0 


3.0 


1-K-Vd 


4.0 


1.0 


1-Q-Vb 


6.0 


3.0 


1-K-VIIIb 


4.0 


1.0 


1-U-VIb 


6.0 


3.5 


1-V-Vb 


5.0 


0.0 


1-R-VIb 


6.0 


4.0 


1-S-Vb 


5.0 


1.0 


lio-H-Vlllb 


3.0 


0.0 


1-C-IIb 


5.0 


2.0 


lio-C-VIb 


5.0 


2.0 


1-C-IIIb 


5.0 


2.0 


lio-T-IIIb 


6.0 


1.5 


1-C-IVb 


5.0 


2.0 


lio-W-VIIIb 


6.0 


2.0 


1-C-Vb 


5.0 


2.0 


25-H-VIIIb 


3.0 


0.0 


1-C-VIb 


5.0 


2.0 


25-K-VIIIb 


4.0 


1.0 


1-C-VII 


5.0 


2.0 


25-C-Vlb 


5.0 


2.0 


1-C-VIIIb 


5.0 


2.0 


25-C-VIIIb 


5.0 


2.0 


1-T-IIIb 


6.0 


1.5 


25-T-IIIb 


6.0 


1.5 


1-T-Vb 


6.0 


1.5 


25io-H-VIllb 


3.0 


0.0 
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^Arabic numbers refer to the type of chemical evolution model (sec text and Tosi 1988a). 
Capital letters indicate the adopted primordial abundances. Roman numerals indicate the stellar 
yields adopted for ^He: I, as described in section 2; II, as for Model I but with g3=l for M < 2 Mq; 
III, as for Model 1 but with g3=l for stars with 1<M/M0 <2; IV, as for Model 1 but with total 
^He destruction in stars with M < 1 Mq; V, following Pinsonneault (1995); VI, following Hogan 
(1995); VII, as for Model VI but with ^He destruction for M < 1.6 Mq; VIII, as for Model VI but 
with ^Hc destruction in the range 1.3 <M/Mq <2.5. Lower case letters indicate the adopted infall 
abundances: a, X2i„/=X2p and X3j„/=X3p; b, X2i„/=0.8X2p and X3j„/=X3p; c, X2i„/=0.8X2p 
and X3j„/=1.2X3p; d, X2i„/=0.8X2p and X3j„/=0.8X3p; e, X2i„/=0.7X2p and X3,„/=X3p. 
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to = 13 Gyr = 10 Gyr 
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to = 13 Gyr = 10 Gyr 
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